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Coating a three-dimensional topological insulator with a ferromagnetic film supporting an in-plane 
magnetic vortex renders electronic states bound to the vortex. Charged zero-modes also appear as 
fully polarized current localized near the vortex center. From the magnetic point of view, such a 
special current can be faced as an alternative technique for detecting magnetic vortices in thin films. 
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I. INTRODUCTION 

Topological Insulator (TI) is a nontrivial class of band 
insulator characterized by a topological index associated 
with the Bloch wave-function in the momentum space. 
They are materials with a bulk band gap generated by 
strong spin-orbit coupling and topologically protected 
metallic surface states. At present, such a behavior 
has been observjed in Bi^Shi-x alloys, Bi2Se^, Bi2Te^, 
Sb2Se3 crya.talsErl3 and in ternary intermetallic Heusler 
compoundsQiJ. The reason for the considerable interest 
in topological insulators lies in their potential for pro- 
ducing new physical phenomena as well as future tech- 
nological applications. These rely on the special prop- 
erties of the topologically protected electronic surface 
states, composed of an odd number of Dirac fermions, 
familiar to the physics of grapheneQ. The physical man- 
ifestation of the topological order comes in the form of 
protected, gapless surface states that are robust against 
damaging the surface by chemical or mechanical means 
such as alterations in its shape or orientation with re- 
spect to the crystal lattice, or even by disordering tl, 
bulk, as long as such changes are applied in moderatior 
The above statements apply exclusively to systems that 
respect time- reversal-symmetry (TRS), that is, nonmag- 
netic TI in zero external magnetic field. When TRS is 
broken, even by a weak magnetic perturbation, a gap can 
be open up into the spectrum of the topologically pro- 
tected surface states. When a gap opens as a consequence 
of a magnetic perturbation, the resulting surface is not an 
ordinary insulator; instead, it is a quantum Hall insulator 
with properties similar to those of the familiar quantum 
Hall systems realized in two-dimensional electron gasS. 
A tunable energy gap at the surface Dirac point provides 
a means to control the surface electric transport, which 
is of great importance for applications. 

Topological order in proximity to magnetism has 
considered as one of the main interest in the fieldnH 



In principle, it can be realized by several means such as 
the application of weak magnetic fields through a Zee- 
man termtfl Hz = —gfiBS ■ through exchange cou- 
pling to magnetic i,ljm _BJmsliil Hint = — AAf • s, and 



magnetic impuritieg3ll3't3 Himp — ' sS{f — Ri). 

The exchange coupling arises due to the spin-spin in- 
teraction between spins in different atoms. Such an in- 
teraction is due to the superposition (overlap) between 
the wavefunction that describes conduction electrons on 
the surface of a TI and bound electrons in the mag- 
netic impurity on the surface. This is one of the most 
appealing from the theoretical point of view because it 
not only gives a simple mechanism to develop the the- 
ory of the topologically quantized magnetoelectric term, 
but it also allows us to look for unexppeted effects that 
can alter the magnetization dynamicalailj. Therefore, 
here we consider a 3D TI coated by thin film of layered 
two-dimensional {2D) Heisenberg ferromagnet. Indeed, 
magnetic materials (ferromagnetic or even antiferromag- 
netic) in two spatial dimensions may support topologi- 
cal excitations such as skyrmions and vortices. Partic- 
ularly, vortices arise in classical magnetic systems con- 
taining an easy-plane anisotropy, which makes the spins 
to prefer to point along the XF-plane. For instance, 
easy-plane ferromagnets are described by the Hamilto- 
nian H = -JY.^,j[SfS'^ + Sf^J + A^f^J], where J > 
is the exchange constant, < A < 1 is the easy-plane 
anisotropy and Si = {Sf, Sf , Si) is the classical spin vec- 
tor at site i. The spin field can be parametrized by two 
scalar fields cj) and e — cos 6 (the azimuthal and polar 
angles) as follows: S = — cos(/),-\/l — sm(j),e}. 
Taken into account the most realistic cases of discrete 
lattices, and depending on the range of A, such an easy- 
plane system supports two types of static vortices: the 
in-plane V|OJjtex (in which all spins are confined to the 
Ary-planetZl) and the out-of-plane vortex (in which some 
spins around the vortex center can point perpendicularly 
to the XF-plane) . Indeed, considering a critical value of 
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A denoted by Ac, then, for the range A < Ac, the stable 
excitation is the in-plane vortex wlaije, for A > Ac, the 
out-of-plane vortex becomes stablet3. The stability of 
these solutions has only been determined via computer 
simulations. The critical anisotropy Ac depends on the 
lattice geometry: for the square lattice, Ac = 0.72; sim- 
ilarly Ac = 0.86 for the hexagonal lattice and Ac = 0.62 
for the triangular latticell3. Qualitatively, similar results 
can be obtained for 2D easy-plane antiferromagnetic sys- 
tems. 

Our interest here is to investigate how an in-plane 
ferromagnetic vortex affects the electronic states lying on 
the surface of a 3D TI. Although 2D magnetic materials 
can support vortex excitations, there is, however, the 
experimental challenging of getting aiid_confining a 
single static vortex in these thin filmsEjEil. Even the 
observation of these excitations in 2D layered easy-plane 
magnetic materials is not a simple task: Indeed, only 
indir aci j. evid ences for their appearance has been reported 
so farll2l'E3l23. In addition, there is the problem of getting 
a static single vortex living in the system for a long 
timeEil. To overcome these difficulties, a possibility 
should be the use of a vortex pinning mechanism. In 
2D magnetic materials, a nonmagnetic, impurity (a 
vacancy) is able to pin tbfi.[jap±eK centeroo, localizing 
the excitation around itE3EZIE3'cj. Another possibility 
(perhaps the simplest) should be through the use of 
thin magnetic nanodisks over the TI. Really, depending 
on the size of a nanodisk, magnetic vortices are nat- 
urally the ground state of these systems and thejf-tjajri 
be directly observed by experimental techniquesEllES. 
Moreover, the vortex ground state is always out-of-plane 
(in the vortex core, the spins point perpendicularly to 
the disk plane) due to the competition between exchange 
and dipolar interactions. However, the presence of 
artificial defects pintantionally inserted in the nanodisk 
(such as holes)E3c3, may destroy the usual out-of 
plane spin structure in the vortex core (after vortex 
capture), resulting in anj-cssentially confined planar 
vortex centered at the hola23. 

Here we are interested in the in-plane vortex effects on 
the TI. Then, in principle, it can be investigated either 
by applying a thin coating of magnetic material (with 
A < Ac and containing a small percentage of nonmag- 
netic impurities) on the TI or by putting a magnetic 
nanodisk (with a hole, preferentially in the disk center) 
on the top of the TI. In the context of the first situa- 
tion, recently, the authors of Ref. [|l3| have theoretically 
studied the magnetic textures, such as domain walls and 
vortices in a ferromagnetic thin film deposited on the sur- 
face of a TI. These authors have found that the magnetic 
textures interacting with the topological surface states 
are electrically charged. The experimental observation 
of massive Dirac Fermion on the surface of a magneti- 
cally doped TI was reported in Ref. ||] , where the authors 
introduced magnetic dopants into three-dimensional TI 
8128 and observed the massive Dirac fermion state by 



angle-resolved photoemission. Other aspects as Hedge- 
hog spin textuseEj, magnetic scattpiug of the carriers 
on the surfaceEa, magnetic,, orderingr'^r^l and topological 
magneto-electric effectEElEilcj were also investigated. 

Considering an in-plane closed fiux magnetized state 
interacting with a 3D TI, we have realized that the inter- 
action between the surface carriers and the vortex leads 
to bound-states and spin polarized zero-modes, besides 
scattered states. Such a framework provides a mechanism 
to obtain spin polarized current based on the zero-mode 
states; conversely, such a physical mechanism appears to 
be useful for experimental magnetism to precise and effi- 
ciently detection of magnetic vortices in layered magnetic 
thin films. Our article is outline as follows: in Section II 
we present the model which takes in to account the basic 
properties of the TI surface electronic states along with 
their coupling to the magnetic moments (spins) lying on 
the ferromagnetic layer. Section HI is devoted to the 
zero-mode states induced by the vortex. Besides its ba- 
sic features we also point out possible utilities for such a 
special current. In Section IV the remaining bound-states 
are presented and discussed, completing the 'atomic pic- 
ture' provided by the compound system vortex- TI charge 
carries. We finally close our paper by pointing out our 
conclusions and prospects for future investigation. 

II. THE MODEL 

The low energy excitations on the surface of a TI are 
described by a Dirac-like equation in a two-dimensional 
plane. When the chemical potential is confined between 
valence and conduction bands (in the ground state), the 
positive energy modes (that describe free conduction 
electrons) are empty, while the negative energy modes 
(that describe bound valence electrons) are filled with 
valence electrons. The mass term in the equation 
expresses the gap of energy between the conduction 
and valences bands (a trivial insulator) and a vanishing 
mass term leads to., a gapless spectrum conductor 
or TI surface modesCJ. The latter ones offer a unique 
platform to investigate the physics of robust Dirac points. 

Here, we analyze the simplest case of a single, Dirac 
point (for instance, Bi2Se^, Bi2Te^ and Sb2Te^S). Let 
us consider the surface of the 3D TI in the ATF-plane (its 
bulk extends to 2: < region) coated by a ferromagnetic 
thin film supporting an in-plane vortex. Therefore, the 
low-energy effective Hamiltonian for the TI surface states 
coupled to the vortex magnetic moments (spins) reads: 

H = Ho+ H,nt (1) 

where 

Ho = i-F^t(r)(fc,a^ -I- kyayWr) , (2) 

and 

i?mt = - ^^M^{r^■m■ (3) 
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Above, s{r) = tp^r)a'ip{r) is the spin density of sur- 
face electrons in the position r, M(r) is the magnetic 
moment apphed on the surface of the TI, which charac- 
terizes the magnetic configuration on the surface and A.; 
is the ferromagnetic exchange constanilJ^B, accounting for 
the coupHng between the electrons lying on the TI surface 
and the magnetic moments composing the magnetic ma- 
terial (exchange- type proximity effect; note that Hint is 
minimized when the electrons spins align to the dipoles) . 
For instance, regarding the material Sb2Te3, doped with 
vanadium the surface exchange parameters A^,, A^, A^ 
are estimated to be of the order of 0.1 — jQ.Sel/, relying on 
the overlap between the wave- functional. For the case of 
a ferromagnetic ordered state perpendicular to the sur- 
face, the interaction Hamiltonian (|^) has the form of a 
mass term for the surface electrons proportional to cTz, 
being the carriers massive with a dispersion relation like: 



E = ±vpJkl + kl + {imvfY , 



(4) 



where the positive and negative signs describe the con- 
duction electrons and the valence electrons (or holes), 
respectively. Whenever this is in order a gap opens. 
Eg — 2^v'j, rendering TI to be a quantum Hall insu- 
lator. As we shall show below, the magnetic vortex also 
breaks TRS in the surface, but does not induces a mass 
gap, like above. As consequences we shall find electronic 
bound-states and zero-energy modes around the vortex. 
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FIG. 1. (Color online) Examples of ferromagnetic vortex-like 
patterns with unity charge, g = -1-1, and with profiles (po = 
+tt/2 (curly- type) and (po — (radial-like). The depicted 
arrows represent the magnetization vector field, M. 



Thus, the eigenvalue problem Hip — Etp can be solved 
analytically by the usual ansatz: 



-!/;(r) 



m = 0, ±1, ±2, 



(8) 



A. Magnetic vortex configuration 

An in-plane vortex in the continuum limit may be de- 
scribed by: 



0, 0: 



q arctan 



y-yo 



(5) 



where q accounts for its topological charge (vorticity), 
{xo,yo) is its centered position while 00 is a constant 
determining its profile (in Fig.^ we depict two exam- 
ples of unity charge vortex, q = I, with distinct pro- 
files). For such a planar vortex, its magnetization, A'l, 
is essentially (in polar coordinates): M{r) — Alcf) with 
= + My — constant and therefore, the interac- 
tion Hamiltonian (Bl) becomes: 



where the components of the spinor ip{r) and xi''') satisfy 
the first order coupled differential equations given by: 



d m+1 M \ , , . E . . 

:r + 7— = 'i—'pi^) ' (9) 

dr r nvp J nvp 



d m M \ , ^ E , ^ , , 



Such expressions may be reduced to a second order 
decoupled differential equation, but to analyze zero 
modes solutions, this form is more appropriated, as 
below. 



III. ZERO-ENERGY MODES 



H,nt^-^\r^M$-dijir}, (6) 

with M — AiMi (hereafter we normalize A^ = IVi). 
This term does not open a mass gap in the spectrum, 
like an ordered ferromagnet does; It rather induces the 
appearance of bound-states, including zero-modes, in the 
electronic spectrum of the composed system, as we shall 
see below. Thus, the dynamics of the low energy states 
is described, in polar coordinates, by the Hamiltonian: 



Firstly, we investigate the zero energy solutions (zero- 
modes) of the equations (||) and ( |lo| ) for the case of a 
planar vortex with 00 = 7r/2. When the Dirac equation 
has a homogeneous position independent mass term, the 
solutions are as usual: continuum solutions with posi- 
tive {E > 0) and negative {E < 0) energies. However, 
in the presence of a topologically non-trivial background 
or defect, a mass term position dependent can appear 
and the spectrum contains (in addition to the habitual 
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positive and negative energy modes) also isolated zero- 
energy mode states. The zero-energy solutions of the 
Dirac equation do not require a specific form for the 
inhomogeneous mass profile; instead, they only belong 
to a nontrivial topological class. Then, the presence of 
zero-energy modes can be established, a priori, by math- 
ematical index theorems, which relate the occurrence of 
these zero modes to the geometry and topology of the 
space on which the Dirac equation is statedcS. The in- 
dex theorem offers an analytic tool that relates the zero 
modes of elliptic operators (like the Dirac Hamiltonian in 
(2-fl) dimensions) with the geometry of the manifold on 
which these operators are definedcj. This theorem can 
be employed to obtain information about the low energy 
behavior of the TI and, in particular, about its conduc- 
tivity properties. A general studji_of the index theorem 
in TI will be presented elsewherecj. 

Now, we investigate the possibility of occurrence of 
zero-modes on a surface of the TI without a mass profile 
in the Dirac equation due to the presence of a non-trivial 
background for the surface carriers. We study this occur- 
rence based on analytical procedures by explicitly solving 
the Dirac equation. In the case of zero energy solutions, 
the radial components of the spinors and ( p^ ) must 
satisfy: 



+ ^(r) =0 

dr r nvp 



d_ 

dr 



-T—)x{r)=0 



(11) 



(12) 



There are two possibilities, depending on the orientation 
of the local magnetization: When the magnetization is 
oriented clockwise, M < 0, the physical zero-modes read: 



Physically, it is precisely the coupling between the or- 
bital angular momentum of the carriers and the mag- 
netic moments of the vortex that yields polarized charged 
zero-modes confined to a circular region with radius 
rg = hvp l\M\^ around the vortex center. Taking Sh-iTe-i 
as an example, one has hvp ~ 3.7eVA while \M\ is 
typically ^ 50meV giving us rg ~ 74A (see Fig.||). It 
is important to remark that, in this case, there is only 
one fermion per quantum state (which is denoted by the 
quantum number to); indeed, as we have seen, the vortex 
chirality requires that these states are polarized. There- 
fore, in principle, such zero modes may contain a very 
large number of polarized carriers (one per each possible 
to) for a given vortex chirality. 





M < 0, S^, = -hi'l 



M > n, S'„ = +r)/2 



FIG. 2. (Color online) By breaking time reversal symmetry 
a ferromagnetic vortex forces the electrons spins, confined to 
the topological insulator surface, to alig n along its magneti- 
zation direction, according to Eqs. ( |l3| ) and (llj) . Such an 
alignment is fully accomplished by zero-mode fermions that 
appear concentrated in a region of radius ro around the vor- 
tex center, what yields the possibility of highly spin-polarized 
zero-mode current. 



(13) 



where a„i is the normalization constant give by a„ = 

m+l/2|+l/2 

1 ' ^ . Solution (111) can be 



^2TT(\2m+l\\) \ 

normalized only if 2to +1 < 0; in addition, these so- 
lutions are spin-polarized with = — ^ , so that the vor- 
tex provided a mechanism for supporting spin-polarized 
zero-energy current on the TI surface (see Fig. ^) . On the 
other hand, if the local magnetization is oriented coun- 
terclockwise, M > 0, the zero- mode solutions are spin 
polarized with = +| and given by: 



i/;(f) 



(14) 



1^ ) . This solution has 



m+l 



with bm = , 

^27r(2m-(-l!) 

physical meaning only if 2to -t- 1 > 0. Both sets of spin- 
polarized zero-modes carry integer electronic charge: 



Q/e 



= / d^rf = / dVV'V = 1 



(15) 



IV. BOUND STATES 

In addition to the zero-modes at there also exist ad- 
ditional bound-state solutions of Eqs. (^ and ( |To| ) with 
-B < 0. The solutions correspond to internal modes in 
which the electrons wave- functions undergo harmonically 
oscillations around the vortex center. Bound-states can 
be induced in 2D or 3D TPs by vacancies or holes in the 
band gap. These defects can affect the transport proper- 
ties of TPs via additional bound-states near the system 
boundary due to the coexistence ofr-the in-gap bound- 
states near the edge or surface statescZl. Here, this effect 
arises due to the exchange coupling between the spin of 
the surface carriers and the magnetic moment of the vor- 
tex in the surface. 

To obtain the bound-states we proceed in the usual 
way by looking for solutions of Eqs. (^ and (^0|) with 
negative energy, E < Q. Decouple these equations we 
obtain the following second order differential equations 
for the up component of the spinor (down component 
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satisfies a similar eq., as below): 



(Pip{r) 1 d(p{r) 



dr 



ip{r) =0. (16) 



To solve this equation, we use the correspondence be- 
tween this problem and the 2D quantum Kepler problem 
(2D hydrogen atom). Such an equation admits solutions 
for electrons restricted to move in a plane around the nu- 
cleus, due to the potential 1/r witha = x'^ + y'^. The 
energy eigenvalues for the bound-states read: 



E^{n, m) = 



(m + 1/2)^ 
{n+\m\ + 1/2)2 



(17) 



with n = 1 , 2 , 3 , . . . and |m| = , 1 , 2 , . . . n - 1. The 
radial eigenfunction cpi be obtained in terms of the hy- 
pergeometric functionua: 



Here, a„ - 



(2|m|)! 



(n -I- |m| - 1)! 



(2n — l)(n — |m 

(a„.™r)l™le-""'-/2 x 
iFi(-n+ |m| + 1, 2|m| ^ 

M(2m+1 



D! 



1/2 



1, an,mr) ■ (18) 
and iFi is a confluent hyperge- 



ometric function. The component Xn,m{T) of the spinor 
can be obtained from the </?n,m('') by using equation (||). 



E>0 

Cfrcc/scattcrcd states) 



E=0 (jero-modes) 



gap 



E-- 0.9428 M 

(bound- states) 



E= — M (jroiind- state) 



FIG. 3. (Color online) The 'atomic' energy levels allowed 
for zero-modes and bound-states in the composed vortex- 
topological insulator system. An almost continuum spectrum 
is allowed between — |M| and « — 0.94|M|, where the most 
energetic bound-state occurs, while 5 ~ +0.94|Afj gaps this 
state from zero-modes. 



Note that the infinite number of energy levels goes from 
|M| (the ground state) to -2v^|Af|/3 « -0.9428|M| 



(the most excited state) but, in contrast to the 2D Hydro- 
gen atom, the ground state is characterized hy n — go 
while the most excited state is characterized by n = 1 
(Fig. H illustrate the energy levels for such an 'atom'). 
The Bohr radius is ap = ^itj^; which is on the order of 

148 A for typical values of and |Af |. The most proba- 
ble radius for these bound-states are (n -I- l/2)ao (let us 
recall that oq is much smaller than the size of typical mag- 
netic nanodisks containing a vortex ground state. This 
is particularly important if these nanoscaled systems are 
employed for TI coating). In addition, the most excited 
bound-state radial wavefunction are: 



Xioir) = ± 



3i / 2M m\ 



\M\r/hvF 



(19) 



(20) 



from what becomes more evident that the effective 
potential brought about by the magnetic vortex is able 
to confine the TI surface electrons in a limited region 
around the vortex center. Similarly to the zero energy 
modes, the region of large probability is a circle with a 
radius on the order of hvF/\M\, depending only on the 
value of the local magnetization M, decreasing as \M\ 
gets higher. In addition, bound-states with normalized 
wavefunctions are possible only if M{2m + 1) > 0. 
Therefore, there are two possibilities: if Af > then 
m > —1/2, while for A/ < 0, m < —1/2 must be taken. 
As a consequence, all bound electrons move around the 
same sense, determined by the vortex magnetization 
chirality (see Fig. ||). 

Besides the bound-states discussed in the former sec- 
tions, a continuum spectrum of scattered states should 
takes place; they could be parametrized by a continuum 
wave-vector k whose associates wave-functions being so- 
lutions of equation (|l^). For these states, the relevant 
results can be obtained in terms of the phase-shifts of 
the carriers, crosa^ection and resistivity, and will be pre- 
sented elsewhereCj. 

Some additional remarks are now in order. In gen- 
eral, the Dirac equation has negative and positive en- 
ergy solutions. The negative energy solutions correspond 
to the states in the valence band while the positive en- 
ergy solutions correspond to the states in the conduc- 
tion band (considering the chemical potential located at 
E = 0). In the ground-state, all the negative energy lev- 
els are filled and, therefore, the ground-state is degener- 
ate (since the zero modes energy are degenerate). Then, 
what is the ground-state at TI surface state in the pres- 
ence of a magnetic vortex? In other words, how could 
one treat the zero-energy modes whenever constructing 
the ground-state of the TI- ferromagnetic composed mate- 
rial? Should the states be empty like the case of positive 
energy or should they be filled as they are for negative 
energy states? 
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If they are filled (or partially filled) when an elec- 
tric current is established over the TI with a vortex 
(and only zero- modes are excited), then the current 
is spin polarized, being the polarization given by 
sense of the magnetization. However if the bound- 
states are excited to the conduction band (overcoming 
S = 2V2|M|/3 « 0.9428|Af|), then the current has a gap 
and the current will display steps when increasing the 
external electrical field similar to the Hall effect under 
increasing of the external magnetic field. 

If the previous study were performed for a planar 
radial vortex with magnetization exhibiting radial-like 
profile, (/)o = 0, we should have realized that the zero- 
modes were not spin polarized while all bound-states, 
with E < 0, were ruled out. 

It is noteworthy to mention that the authors of Ref. 
pO| have treated a similar problem of a magnetic vor- 
tex coating the surface of a TI. However, in their work 
they have considered the question of whether an effective 
gauge field would be induced by strain on the surface of 
a topological insulator. They have shown that the gauge 
fields arise when the surface is coated with an easy-plane 
ferromagnet with magnetization parallel to the surfaceE^I. 

V. DISCUSSIONS AND CONCLUSIONS 

In summary, we have investigated the effects of 
a closed fiux magnetization (in-plane vortex with 
•^0 = ""/S) on the surface of a 3D strong topological 
insulator with one Dirac cone. We have shown that 
the interplay between topologically protected surface 
states and the broken of time-reversal symmetry by an 
in-plane magnetic vortex coating TI surface does not 
generate a mass gap; It rather leads bound-states and 
localized spin polarized charged zero-modes. Indeed, 
considering the bound-states, the very presence of the 
vortex generates an effective potential for the TI surface 
carriers which is equivalent to the problem of a two- 
dimensional hydrogen atom. The range of energy levels 
of the bound-states is very short and goes from — | M | 
to — 0.9428|M| (comprising infinite possible states in 



this range), which is energetically separated from the 
zero-modes, = 0, by a large gap S sa 0.9428|M|; 
a picture of such energetics may be seen in Fig. |^). 
Zero-modes and bound-states decay exponentially as 
exp(~\M\r/hvF) away from the vortex center, which 
means that some electrons are essentially localized in a 
circular region with radius tq = Hvp/\M\ (around 70 A 
for typical TPs with "atomic" binding energy ~ 50meV. 
Physical systems supporting stable zero-modes deserve 
special attention, for instance, in topological quantum 
computation for the ground-state degeneracy remains 
unaffected by small perturbations, what can render error 
free quantum information encoding. 

Instead of large magnetic films, nanoscale magnets can 
also be used to form the composed system with TPs. In 
this case, special care must be taken because vortices in 
confined structures must have out-of-plane components 
near its center, what opens up a gap in the TI electronic 
states. The remaining effects of a vortex, as studied here, 
should remain qualitatively unaffected, still supporting 
zero-modes and bound-states on the TI surface. Even 
in the case of nanomagnets, in-plane vortices may take 
place, for instance, if-arCavity is intentionally inserted at 
the center of the diskEaEil or we employ other geometries, 
such as a nanotorusEll. Another interesting possibility is 
to coat TI surface with an elliptic thin disc supporting a 
pair vortex-antivortex, with opposite chiralitics. On the 
TI surface it would be induced two distinct zero-modes 
fermionic current flowing in opposite directions, around 
each of the vortex center. 

Finally, it should be remarked that the detection 
of the electronic modes predicted here on the surface 
of a magnetically coated TI could be an interesting 
alternative method to experimentally detect vortices in 
layered 2D Hcisenberg ferromagnets. 
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